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Abstract 

The optimal conformal mapping of the Borel plane was recently used to acceler- 
ate the convergence of the perturbation expansions in QCD. In this work we discuss 
the relevance of the method for the calculation of the Laplace-Borel integral ex- 
pressing formally the QCD Green functions. We define an optimal expansion of the 
Laplace-Borel integral in the principal value prescription and establish conditions 
under which the expansion is convergent. 



1 Introduction 



A suitable method for accelerating the convergence of power series is based on conformal 
mappings. As is known, a power series converges inside the circle passing through the 
nearest singularity of the function to be approximated. Some time ago, in Ref. fl|], it 
was shown that, if the position of the singularities of the expanded function is known, 
one can reach the fastest convergence rate by expanding in powers of the function that 
conformally maps the whole holomorphy domain onto a unit disk. In addition, the conver- 
gence region extends over the whole holomorphy domain. In a recent paper 0] we applied 
the technique proposed in [JTJ] to the Borel transform of the Green functions in perturba- 
tive QCD. As discussed in 0, the Borel plane is very suitable for applying the method, 
since some information about the singularities of the Borel transform is available from the 
study of certain classes of Feynman digrams and from nonperturbative arguments. By 
the technique of conformal mapping, this additional information can to a certain extent 
be incorporated even into the lowest-order terms. In this way, the convergence of the 
perturbative expansion is improved, allowing one in particular to approximately predict 
the next-order perturbative terms from the calculated low-order ones ||. 

In Ref.0 the expansion in powers of the optimal conformal mapping variable^] was 
also used to calculate the Borel-Laplace integral, which is supposed to give, with a certain 
prescription of treating the infrared renormalons, the Borel summation of the large orders 
in the Green functions. The numerical results on mathematical models discussed in || 
indicate that the power expansion in the optimal variable makes also the calculation of 
the Borel integral convergent, in addition with a very high convergence rate. However, 
only qualitative arguments explaining the results were given, and the problem whether the 
improved expansion of the integral is convergent, or signs of divergence might appear at 
large orders, remained open. In the present paper we address this problem and investigate 
the convergence of the expansion of the Borel integral in perturbative QCD, improved by 
the use of conformal mapping. 

2 Optimal expansion of the Laplace-Borel integral 

We study the following integral 

oo 

1(a) = J e-«B(u)du, (1) 
o 

where B(u) is assumed to be analytic near u = 0, where it can be expanded as a Taylor 
series 

oo 

B(u) = Y,b n u n (2) 

n=0 

lr rhe conformal mapping that maps the whole holomorphy domain of the expanded function onto the 
unit disk will be called optimal. In this case, the singularities are mapped onto the boundary circle, and 
the requirement of holomorphy implies convergence of the power series at every point of the disk, which 
is the map of the holomorphy domain. 
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converging inside a circle of non- vanishing radius. The function 1(a) is of interest for the 
Borel summation of the Green functions in perturbative QCD; note, however, that the in- 
tegral on the right hand side of (|l|) is ill-defined if B(u) has singularities along the positive 
real semiaxis, which is the case of QCD (infrared renormalons) . First, the singularities 
of B(u) (renormalons of either kind) make the expansion (|2]) badly divergent along the 
integration path. Second, the function B(u) itself is, because of infrared renormalons, not 
uniquely defined along the integration path. We shall discuss both these problems below 
in this paper. 

We shall consider for illustration the Adler function D(s) of the massless QCD vacuum 
polarization, which can be expressed formally as @-|| 

D(s) = 1 + -J- 1(a) , (3) 

with a = Poa s (—s), where a s (— s) is the running coupling, and (3q = (33 — 2n/)/127r 
is the first coefficient of the (3 function. The expression (Q) formally reproduces the 
renormalization-group-improved expansion of the Adler function 

71 



^I^ + E^Pt 2 > ( 4 ) 



71=1 

by taking the coefficients b n in the expansion (0) of the form 

1 Dn+l 



(5) 



We consider also minkowskian quantities, like the hadronic decay rate of the r lepton, 
R T , which can be expressed formally as || 



R T = 3(1 + 5 EW ) 



1 + JL r due xp ( - U ) B(u) F(u) 
tt/3 Jo \ f3oa s [mi)J 



(6) 



Here Sew is an electroweak correction, B(u) is the Borel transform of the Adler function 
and 

F(u) = ~ 12 sin(?rM) (7) 

The extra factor sin(7TM) in the Laplace-Borel integral is generic for minkowskian quan- 
tities. Note that strictly speaking the expressions @ and (|]) are not equivalent to the 
Borel summation method, which requires an analytic continuation of B(u) from the con- 
vergence disk to an infinite strip of non-vanishing width, bisected by the real positive 
semiaxis. This condition is not fulfilled in QCD because of infrared renormalons, which 
produce cuts of B(u) located along the real positive semiaxis. 

We shall be concerned with the evaluation of the integral ([I]) for complex a of the 
general form a = |a|e 1 ^, where ip = arg a is the phase of a. In the case of the Adler function, 
with the running coupling at one loop in the V scheme a^^(— s) = l/[/?o hi(— s/Ay)], and 
writing — s = \s \ e*^ _7r \ we have 

I=lnM_ i ( 7r _0). (8) 

a l\y 
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Outside the Landau region, i.e. for \s\ > Ay, we have cosip > 0, so that 



<l (9) 



and i[) is related to the momentum plane variable s by 



ip = Arctg 



[n - 0)7^-2- 



A 2 



(10) 



The phase ip is positive for s in the upper half of the s- plane, where < <fi < it, negative 
for s in the lower half-plane, where n < <p < 2tt, and along the euclidean axis. 

For the minkowskian quantity (|6|) we combine the additional factors exp(±i7rw) due 
to the sinus with the exponential, which amounts to taking a complex with 

V> = ± Arctg [ira] , (11) 

where a = p a s (m 2 ). 

As already mentioned, the Borel transform B{u) has singularities in the complex plane, 
correlated to the factorial increase of the perturbative coefficients of the Green functions 
at large orders 0, M. The precise form of the singularities is not known for the exact 
theory, but the position and the nature of the first renormalons can be inferred from 
general principles. In the case of the Adler function the first ultraviolet (UV) renormalon 
is situated at u = — 1 and the first infrared (IR) one at u = 2, and they are branch points 
of the type (1 + w) 71 and (2 — u) 72 respectively, with 71 computed in Ref. pi and 72 in 



J/J. We mention also that the summation of the one- renormalon chains in massless QCD 
in the large (3q limit gives [§], 0: 

32e~ CM ~ (-l) k k . . 

B{u) = — — r y i — ? , 12 

1 ; 3(2 -u) t 2 [* 2 - (1 - «) 2 ] 2 ' 1 ^ 

i.e. all the singularities are poles (C is a scheme-dependent constant, with C = —5/3 in 
the MS scheme, and C = in the V scheme) 0, @. 

The series (@) converges only inside the circle \u\ < R passing through the nearest 
singularity (R = 1 for the Adler function). Since the integration range in ([I]) extends far 
outside this region, by inserting (fj) in (fj) and integrating term by term one obtains a 
divergent expansion. By the technique of conformal mappings one extends the domain of 
convergence of a series beyond the limit imposed by the first singularity. In [0 we used 
the optimal conformal mapping 



y/T+U- Jl-u/2 

w = w(u) = — — = , (13) 

VTTu + Jl - u/2 



with the inverse u(w) = 8w/(3 — 2w + 3w 2 ). The transformation ( p!3|) preserves the origin 
and maps the complex u plane, cut along the real axis for u > 2 and for u < — 1, onto the 
interior of the circle \w\ < 1, all the singularities of the Borel transform being mapped 
onto the boundary \w\ = 1. The expansion in powers of w , 

00 

B(u) = Y / c n w n , (14) 

n=0 
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is called optimal because it converges inside the circle \w\ < 1, i.e. in the whole domain of 
holomorphy of B(u) (which is the doubly cut complex w-plane in our case), up to points 
close to the branch cuts produced by renormalons. As was already pointed out above, this 
power expansion yields, when compared with other conformal mappings, the fastest large- 
order convergence rate (see a proof in [|I|]). In practice, as discussed in 0, the expansion 
(|i~4]) is obtained by suitably reorganizing the summation of the original series (§). More 
precisely, consider the expansion of each u n in powers of w, truncated at a finite order N. 
In particular, in our case this expansion has the general form 



N 



U 



N 



Y^C n jW 3 , (15) 



with the coefficients c n j obtained by expressing u in terms of w (in our case u(w) is 
given explicitly after formula (|T3|)). Starting now with the expansion (fj) truncated at 
finite order N, and replacing each u n by its approximant u' N , one obtains a truncated 
expansion of the function B in powers of w, which in the limit N —>■ oo gives (fH|). 

By inserting the optimal expansion (|14D into the integral ([!]) we obtain the formal 
development 

oo 

I{0) = E C nW , (16) 
n=0 

with 

oo 

I n (a) = f e~^w n du. (17) 



In the present work we shall adopt (|TB|) as the optimal expansion of the Laplace-Borel 
integral. We point out that in the physical case this seems to be a natural definition. 
Indeed, when attempting to make the Borel summation of a perturbation expansion in 
QCD, one starts with a finite sum of the form 



N f u 

^(a) = E & ™ / e ~ du • 

n=0 n 



By replacing here the powers u 11 with the approximations (|15D , we replace In (a) by an 
expansion of the form 

JV 

E C nW , (19) 

n=0 

which in the limit N — > oo leads to (filj). 

Actually, as mentioned above, the conditions for the Borel summation are, because 
of the infrared renormalons, not fulfilled. Therefore, Eq. (O) can be considered as a 
definition of 1(a), provided that (i) the integration path in expressions like ([!]) or (|ITD is 
consistently defined, and (ii) the series fliTf ) is convergent. Let us devote a brief discussion 
to these conditions. 

(i) As concerns the integration contour, let us notice that the expansion (|16|) has not 
a precise mathematical sense with the I n (a) defined by (|T7D, because the integration path 
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runs along the positive real semiaxis, where the w n have cuts. We shall adopt, as in ||, 
the generalized principal value (PV) prescription, defining the I PV (a) as 

C » = \ J e-S («;(«))" dn + \ J e"2 («;(«))» d M (20) 
c+ c_ 

for n = 0, 1, 2, where C + are lines parallel to the real positive axis, slightly above 
(below) it. While the PV prescription does not always give the expected results |14|], in 



QCD it has the advantage that it reproduces, to a larger extent than other choices, the 
momentum plane analyticity properties of the Green functions derived from the general 



principles of field theory. In particular, as discussed in [15] , the Adler function calculated 
with the PV prescription has no unphysical singularities in the region \s\ > A 2 . The 
functions I PV (a), n = 1, 2... are chosen so as to share some of the known properties with 
the unknown I PV (a). This makes them suited for the definition of I PV (a) by means of 
the expansion 



rPV/„\ „ T PV 



E ^fW • (21) 

n=0 



The convergence of the series (21) for complex a is not a priori obvious. Indeed, 



the expansion (14) converges at points \w\ < 1, therefore in the neigbourhood of the 



integration axis, but not necessarily on the boundary. One might therefore expect that 
the boundary singularities could manifest in a dramatic way for very large orders N, 
making the series (|2~1~D divergent, like in the case of the original expansion (0). In || 
we investigated mathematical models with B(u) having a few number of isolated branch 
point singularities, and real values of a. The numerical results confirm that the expansion 
(0) in powers of u gives results which deviate dramatically from the exact value for large 



N, which is typical for a divergent expansion. On the other hand, the improved series (14) 
in powers of the optimal variable led to results improving continuously with increasing 
N, and no signs of divergence appeared even at very high N. In the next Section we shall 
discuss the convergence of the optimal expansion of the Laplace-Borel integral, bringing 
analytic arguments which explain the numerical results obtained in ||. 



3 Convergence of the optimal expansion 

We investigate the expansion ([H]) with the functions I PV (a) defined by means of the PV 
prescription (|20"1) . We consider in our discussion analytic functions B of real type, i.e. 
which satisfy B*(u) = B(u*), where u* is the complex conjugate of u. Therefore, the 
coefficients b n in the expansion (|J), as well as the coefficients c n in the expansion (|14"D are 
real. 

The contribution to (|20|) of integral along the contour C + can be written as 

/» = Je- F ^du, (22) 

c+ 

where 

u 

F n (u) = - -n\nw(u) . (23) 
a 
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We evaluate the integral fl22|) for large n by applying the method of steepest descent JT2 
TM. The saddle points are given by the equation 



w'(u) 1 
w(u) an ' 

which has four solutions, having at large n the form 



(24) 



1 + z 1 — i — 1 + i — 1 — z 

-^ITTV^ -^v^, 2l/4 v^, 2 i /4 v 7 ^- (25) 
Of interest for the evaluation of ( p2|) is the point 

u = 2" 1/4 (1 +i)v^= |w |e* Q (26) 

with 

K|=2 1 /4 V ^ ; « = ^ + ^> ( 27 ) 

which is situated in the first quadrant of the w-plane. Indeed, since the phase ip of the 
parameter a satisfies the condition (|9]), then Reu > and Imu > 0. 
Near the saddle point F n (u) can be expanded as 

F n (u) = F n (u Q ) + \F' r [{u ){u - u Q ) 2 + ... . (28) 

By using the expansion of w(u) for large u in the upper half plane ( w ~ ((1 — iy/2/u), 
where ( = (y/2 + i)/(y/2 — i)), we obtain after a straightfoward calculation 

e -F n (u ) w ^ _ 2 ' y ^ r _ ] n e - 2 ~ 1/4{1+i) \^ « C e" 23/4(1+i) v 7 ? (29) 

and 



(1 + i)y/cm 



F^uo)* 2 ^^ = |^'M|e*, (30) 



IP'V M 23/4 a * ^ 

ln\a\ 



where 

o3/4 

,„J3' 4 2 

Therefore ( p^) becomes 

^ Ce" 2374 ^^ / e"^-* du. (31) 



In order to evaluate the integral we first rotate the contour C + in the trigonometric 
direction in the upper half-plane, until it becomes a line passing through the origin and 
the saddle point Uo- The rotation is possible since B(u) (and therefore also w) has no 
singularities outside the real axis, and the arc of the circle at infinity gives a vanishing 



6 



contribution, as can be easily verified. Along the rotated line u = e ia t, where a is the 
phase of uq defined in (^?p and t is real, so the integral in (^Tj) becomes 



e m / e 



e i(2a+^) ( t _| uo |)2 



dt . (32) 



Since cos(2a + /5) > for ip satisfying the condition (|^), the integration axis lies in the two 
valleys near the saddle point u$. Therefore it can be deformed into the path of steepest 
descent through u , without passing outside the valleys. We take the integral along the 
path going to infinity 

u-u * ^Me^/) (33) 
with real p. The phase of (u — u ) 2 exactly compensates the phase of F"(u ), making 



the exponent of the integrand in ( |31"D real. The integrand can be written as e p and the 
integral done explicitly gives 



W^Ce-^'+Vf j£, (34) 

SP 2 



i.e. up to a constant independent of n 

J+(a)«niCe- 23/4 ( 1+ ^. (35) 

It is important to note that the path of steepest descent must not cross the real axis, 
where B(u) has singularities. From (j33|) this implies —(3/2 = n/8 + 3ifi/4 > which 
writes as 

V>~- (36) 
o 

The evaluation of the integral along the contour C_ in (|20| ) proceeds in a similar way. 
The saddle point of interest is 



u'q = 2~ 1 / 4 (1 - = 2 1 / 4 J\^e- i ^-* ) , (37) 



which satisfies Rewo > and Imuo < for ip in the range given in (H). Instead of (|3T|) we 
I- (a) « (CTe-^^V? / e-^^' (— o) 2 du , (38) 



where /3' = it/A — 3ip/2 . We rotate the contour C_ in the lower half-plane up to a line 
passing through the point u' Q , and then deform it into the steepest descent path. One can 
easily verify that this path does not cross the real axis for —tt/8 + 3^/4 < which writes 

as 

^ < \ ■ (39) 
o 



Collecting the terms we obtain the coefficients I n (a) in the PV prescription fl2"Dp as 

J n (a) « nk n e- 23/4(W) + n\ {p)"**'^^ . (40) 



7 



In order to examine the convergence of the expansion ([H]), we consider the ratio 



CJ T n{a ) s , (41) 

for large n. If the coefficients c n do not grow too rapidly, i.e. 

\c n \ < Ce e " V2 , (42) 

for all e > 0, then the expansion (|2l"l) converges for a complex in the domain 

Re[(l±i)cr 1/2 ] > 0. (43) 

As we already discussed these conditions are equivalent to @. If the coefficients c n behave 
at large n like 

\c n \^e cnl/2 (44) 
for some positive c, then the expansion (|21"D converges in the domain 

Re[(l ± i)a~ 1/2 + c] > , (45) 

while for coefficients c n which grow faster than exp(cn 1,/2 ) the new series ( pi]) is also 
divergent. We mention that such a behaviour is not excluded in general for series of the 
form fll4]) with a radius of convergence equal to 1 ||14j| . 

We recall however that the expression (pL0|) is valid only for ip which satisfy the condi- 
tions © and (§9|), i.e. 

< | , (46) 

which define a sector in the a- complex plane (we recall that ip is the phase of a). This 
inequality is a condition of applicability of the steepest descent method used by us. We 
found therefore that the expansion (21), improved by the optimal conformal mapping of 
the Borel plane, is convergent if the Taylor coefficients c n of the expansion ( |14|) satisfy 
the condition (fi"2j), at least inside the sector ( fIBf ) of the complex plane of a, or, if they 
behave like (f44|), in the smallest of the domains (f45|) and (0). For the Adler function in 
massless QCD, using (|T0|) we write the condition (§|]) in the form 

''-^Tl" 1 ^ (47) 

where <fi is the phase of s and we have \s\ > Ay. For the minkowskian quantities, from 
(pT|) we obtain 

na < -y= , (48) 

which means in particular that for the r- hadronic decay rate the expansion defined as in 
(13) is convergent for a s {m 2 T ) < 4/(9\/3) « 0.257. 

The behaviour of the coefficients c n depends on the singularities of B(w). By the 
conformal mapping (|13|) all the renormalons are situated on the circle \w\ = 1, appearing 



8 



in conjugate pairs since B(u) is of real type. Assuming that all the singularities are poles 
or branch points, c n has the generic form 

c n « - Re J2 r JPj(Pj + l)(Pi + 2)...(Pj + n)e^^fe + ") , (49) 

where exp(±i/3j) denote the position of the renormalon in the w-plane, rj the residue, 
and pj the exponent of the singularity. In Ref. || we investigated simple models with a 
finite number of singularities, and real values of the parameter a, for which the conditions 
of convergence are satisfied. In the physical case, one knows only that for the first UV 
renormalon at = it and p± = 271, and for the first IR renormalon «2 = and P2 = 272. 
In the large f3 Q case, as seen from fll2D , all the singularities are poles, pj in fl49f) is inde- 
pendent of j, and rj are known. In this case the condition of convergence fl42| ) is satisfied. 
Therefore, the optimal expansion on the Laplace-Borel integral, in the PV prescription, 
for the summation of one renormalon chains in the large (3q limit, is convergent, at least 
in the sector of the complex a plane defined by the condition (|46D . 

In conclusion, we investigated the expansion of the Laplace- Borel integral in pertur- 
bative QCD, improved by the analytic continuation of the Borel transform outside the 
perturbative convergence disk (and, simultaneously, by reaching the fastest convergence 
rate) by means of the optimal conformal mapping [fj]. The convergence properties of the 
new expansion depend on the strength of the singularities of the Borel transform, reflected 
in the behaviour of the Taylor coefficients of the expansion fll4]). If the Taylor coefficients 
satisfy the condition (0), the new expansion of the Laplace- Borel integral converges 
in the sector of the complex plane of the coupling a defined by fl4TS|). The conditions are 
satisfied in the case of the resummation of one-loop renormalons in the large- (3q limit. We 
mention that in the region where the series converges the function 1(a) must be analytic. 
For the Adler function in the complex momentum plane this corresponds to the region 
described by Eq. (|47D , where D(s) is analytic. 
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